Conditional linearizability criteria for a system of third-order ordinary 

differential equations 

o 

O ' F. M. Mahomed'^, I. Naeem" and Asghar Qadir 

> 
O 

00 



> 
in 

(N 



O 



X 



"Centre for Differential Equations, Continuum Mecfianics and Applications 
School of Computational and Applied Mathematics 
University of the Witwatersrand 
Wits 2050, South Africa 
E-mail: Fazal.Mahomed@wits 



< 

u 

] ''Centre for Advanced Mathematics and Physics 
^ ' National University of Sciences and Technology 
^ tt^ ■ Campus of the College of Electrical and Mechanical Engineering 
' Peshawar Road, Rawalpindi, Pakistan and 
Department of Mathematical Sciences 
l/-^ King Fahd University of Petroleum and Minerals 
Dhahran 31261, Saudi Arabia 
E-mail: aqadirs@comsats.net.pk 



Abstract 

We provide linearizability criteria for a class of systems of third-order ordinary differential equa- 
tions (ODEs) that is cubically semi-linear in the first derivative, by differentiating a system of 
second-order quadratically semi-linear ODEs and using the original system to replace the second 
derivative. The procedure developed splits into two cases, those where the coefficients are con- 



, stant and those where they are variables. Both cases are discussed and examples given. 



1 Introduction 



Linearization plays an important role in the study of nonlinear differential equations as there is 
no standard method to solve nonlinear ODEs. In 1883, Lie [T] gave the linearizability criteria for 
scalar second-order ordinary differential equations and he constructed the most general form of 
second-order ODEs reducible to linear ones by changing the dependent and independent variables. 
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Lie also developed point transformation groups [2]. To the best of our knowledge Lie did not 
provide linearizability criteria either for second-order systems or for third-order scalar ODEs, 
leave alone systems of third order ODEs. Nor does this aspect appear in the literature. One can 
use the approximation (perturbation method) to convert nonlinear ODEs into linear ones [3J but 
one faces problems in the convergence of the series. A similar problem arises in numerical schemes. 
Contact transformations were used by Chern [H E] to reduce scalar third-order equations into the 
linear forms u'" = and u" + u = 0. Grebot [B| used fibre preserving point transformations for 
mapping third-order ODEs to linear form. A generalization of this work was later performed by 
Neut and Petitot Ibragimov and Meleshko |8j and Maleshko [9J also studied the reduction of 
scalar third-order ODEs to linear forms. Whereas Neut and Petitot and Ibragimov and Meleshko 
used the classical Lie method, Meleshko dealt with a sub-class of those equations in which the 
independent variable did not appear explicitly and hence one could convert to treating the function 
to be determined as the independent variable and and its derivative as the dependent variable, 
thereby reducing the order by one. 

The linearizability criteria for a system of second order quadratically semi-linear ODEs were 
discussed by two of us [lU], using the connection between the symmetries and the isometrics of 
the system of the geodesic equations [TT]. The criteria require that the curvature tensor be zero, 
i.e. the space be flat. Subsequently [13] a projection procedure of Aminova and Aminov pj| was 
followed, using the translational invariance of the geodetic parameter in the geodesic equation. 
A system of two dimensions was used to get a single cubically semi-linear ODE and Lie's results 
on linearization of scalar second order ODEs was also re-derived. Moreover, invariant criteria for 
a system of two cubically semi-linear second-order ODEs to be reducible to a linear system were 
obtained. 

A method was developed [H] The purpose of this paper is to construct the linearizability condi- 
tions for a system of third-order ODEs, cubically semi-linear in the first derivative by differen- 
tiating the system of second order quadratically semi-linear ODEs and then setting the original 
system to replace the second derivative. The procedure employed here uses the derivative of the 
vector of the coefficients in the equation and solves a system of non-homogeneous algebraic equa- 
tions for it. If this vector becomes zero, the procedure of solving the system by inverting a matrix 
cannot be used directly. As such, one has to go back to the system and solve it directly. This 
allows us to construct the linearizability conditions in all cases. Since there is no work on explicit 
linearizability criteria of systems of third order non-linear ODEs, we cannot compare our results 
with other works. However, there is a result for the number of classes of linearizable systems of 
ODEs, which is based on there being 3n arbitrary constants for the system, using the classical 
Lie method. The class of linearizable equations obtained here is not contained in those classes as 
it is non-classical and does not have the required number of arbitrary constants. 
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The outline of the work is as foUows. In the next section we briefly mention the geometrical 
notation used and review the linearizability criteria for second order quadratically semi-linear 
systems. In section 3 we discuss the linearizability criteria for the general case for a system of 
third-order ODEs when at least one of the coefficients is variable. The special case when all 
the coefficients are constant is discussed in the next section. In section 5 we illustrate via some 
examples to verify our criteria. Concluding remarks are given in section 6. 



2 Preliminaries 



We shall use the Einstein summation convention that repeated indices are summed over. The 
position vector, , {j — 1, ...,n) is assumed to be a function of a geodetic parameter s and ' 

represents the derivative relative to s. The metric tensor will be denoted by gij{x"^) and the 
partial derivative relative to by " ,j" . The Christoffel symbols F^-^ are defined in terms of the 
metric tensor by 

r}fe = 2^'\9ji,k + gki,j - gjk,i)- (2.1) 

The geodesic is given by the equation 

x"" + Tl^x^'x"' = 0. (2.2) 

The Riemann tensor is defined by 

r>i Y^i Y^i i "p^ "pm "pi pm /n q\ 

^jkl — ^ jl,k jk,l "I" mk'- jl ml'- jki \^-'^) 

which satisfies 

R]ki = -Rjik^ (2-4) 

R]ki + ^kij + ^ijk = 0> (2-5) 

and 

^]kl;m + ^]lm;k + ^)mk;l — 0" (2-6) 

The Riemann tensor in complete covariant form is 

Rijkl = 9imR]lkJ (2-7) 

and has the additional property 

Rijkl — —Rjikl- (2.8) 



A general system of quadratically semi-linear second order ODEs can be written as 

x^" + -it^x^'x"' + Plx^' + a" = 0, (2.9) 

where 7^^, a" are functions of the n independent variables. If a" = we will call it of 
geodesic type if the coefficients 7^^ can be written as F^^ for some metric tensor gah- 

For example, a two-dimensional system of geodesic type is of the form 

x" = a{x, y)x''^ + 2b{x, y)x'y' + c{x, y)y''^, 
y" = d{x, y)x'^ + 2e(x, y)x'y' + f{x, y)y''^, 

where we have written x^ = x, x"^ = y and 

Tl, = -a, Tl, = -b, Tl, = -c, (2.11) 

Tl, = -d, Tl, = -e, Tl, = -f. (2.12) 

The linearizability condition of [TO] is that the curvature tensor constructed from these coefficients, 
regarded as Christoffel symbols, be zero, i. e. i?*^; = 0. In terms of the coefficients of (3.2) we 
have 

ay-b.^ + be-cd = 0, (2.13) 
by -c, + (ac - b^) + {bf - ce) = 0, (2.14) 
dy - - {ae - bd) - {df - e'^) = 0, (2.15) 
{b + f). = ia + e)y. (2.16) 

3 Linearization for third-order ODEs 

To obtain the criteria for the third order semi-linear system of ODEs we differentiate (2.2) to get 

x""'" + 2rix^"x'' + ri yx^'x'^' = 0. (3.1) 
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This is a total derivative system and may seem to be an artificial way of getting the third order 
systems. However, if we now replace the in (3.1) using (2.2) we get the system 

x""'" + - 2Tl^raiy' x'^' x'' = 0, (3.2) 

which is not a total derivative system. Thus the general form will be a third order system of 
ODEs that is cubically semi-linear in the first derivative 

x"'" + Al^^x^' x^' x'^' = 0, (3.3) 

which will be linearizable if we can make the identification 



In the case of a system of two third-order ODEs the linearizability conditions (2.13) - (2.16) are 
satisfied. Using (3.3) and (3.4) with the help of (2.13) - (2.16), we obtain 

All = -ax - - 2bd, (3.5) 



^112 ~ '^y 2&a; 


-2{3ab + 2be + cd), 


(3.6) 


^122 ~ + Cx) — 


2(ac + 26^ + 6/ + 2ce), 


(3.7) 


^222 — ~^y 


-2{bc + cf), 


(3.8) 


^111 — "a; 


— 2ad — 2de, 


(3.9) 


^112 ~ ~i^y "1" ^e-r) — 


2{2bd + ae + 2e'^ + df), 


(3.10) 


^122 ~ "(2% + fx) 


-2{cd + 2be + 3ef), 


(3.11) 


^222 — ~fy 


- 2ce - 2f. 


(3.12) 
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Invoking (2.13) - (2.16) we can write (3.5) - (3.12) as 

= -(Ajii + 2a^ + 26c/), = -(^^22/3 + 206 + 266), 

= -(^1^2/3 + cd + 2ab + be), by = -(^^22/3 + 6^ + ac + ce + bf), 

= -(A}22/3 + 26^ + 2ce), = -(A^22/3 + 26c + 2c/), 

4 = -{Al-^^ + 2ad + 2de), dy = -{Al^^/3 + 2bd + 2e^), 

= -(.4^^2/3 + df + bd + ae + e^), Cy = -(^^22/3 + cd + be + 2e/), 

fx = -(A?22/3 + 26e + 2e/), /, = -(A222 + 2ce + 2f). (3.13) 

Using the compatibility conditions and for convenience writing 

^111 ~ ^112 ~ Q' ^122 ~ -^5 ^222 ~ 



Al,,=T, Al,, = U, Al,, = V, Al,, = W, (3.14) 
the system ( 13.131) reduces to 

3Py - 2aQ - 2dR = Q^- 6bP - 2eQ, 
Qy-aR- 3dS = - 3cP - fQ, 
Ry - 2bR - 6eS = 35^. - 2cQ - 2fR, 
3Ty - 2aU - 2dV = - 6bT - 2eU, 
Uy-aV - 3dW = V^- 3cT - fU, 
Vy - 2bV - 6eW = 3W^ - 2cU - 2fV. (3.15) 



For n = 2 the system (3.4) is hnearizable if it satisfies conditions (3.13) and (3.15), where 
a, b, c, d, e and / are given in the Appendix. Our procedure puts the derivatives of the A^j/^j^ 
as a vector on the right side of a system of hnear equations and solves for the constraints on the 
coefficients by inverting the matrix of coefficients on the left side. If the vector on the right is 
zero, we cannot use this procedure to obtain the conditions on the coefficients. As such, for this 
procedure to work at least one of the coefficients in the equations must be variable. If all the 
coefficients in system (3.4) are constants then the derivatives of A'^jj^^ are zero. For this case we 
need to develop a separate procedure that is done in the next section. 



4 Special cases 

In the case Aj,^i ^ = 0, the system (3.15) cannot be used for obtaining the coefficients by inverting 
the matrix. We have to solve the equations without resorting to that procedure. This system 
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reduces to 



~aQ-dR + 3bP + eQ = 0, -aR - 3dS + 3cP + fQ = 0, -bR - 3eS + cQ + fR = 0, 
-aU -dV + 3bT + ef/ = 0, -aV - 3dW + 3cT + fU = 0, -bV - 3eW + cU + fV = 0, 

and equations (3.5) - (3.12) become 

Al^-^ + 2{a^ + bd) = 0, A\-^2 + 6{ab + be) = 0, 
Aji2 + 3(crf + 2ab + be) = 0, + 3(6^ + ac + ce + bf) = 0, 

A22 + 6(&' + ce) = 0, ^^22 + 2{bc + of) = 0, 

Al^-^ + 2{ad + de) = 0, Al-^^ + 6(M + e^) = 0, 
^112 + 3(c^/ + bd + ae + e^) = 0, ^^22 + 3{cd + be + 2ef) = 0, 

^^22 + 6(6e + ef) = 0, ^^22 + 2(ce + Z^) = 0. 



(4.1) 



Non-trivial solutions of (4.1) exist if 



-2Q 


6P 





-2R 


2Q 





-R 





3P 


-3S 










-2R 


2Q 





-6S 


2R 


-2U 


6T 





-2V 


2U 





-V 





3T 


-3W 





t/ 





-2^ 


2U 





-6W^ 


2y 



= 0. 



This gives the linearizability criteria of (3.4) when all the coefficients are constants. 



From system (14. ip . we have 



cd — be = 0, b'^ + ce — ac — bf = 0, bd + e^ — df — ae = 0, 
Al^^ = -2{a^ + bd), A\^^ = -6{ab + be), 
AI22 = -6(&' + ce), ^^22 = -2(6c + c/), 
Aj^^ = -2{ad + de), Aj^^ = -Q{bd + e^), 
^^22 = -6{be + ef), Al,, = -2{ce + f'). (4.2) 



In order to solve system fl4.2l) . the following cases are considered. 
Case 1: 6 = 0. 

The following subcases arise. 
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Case 1.1: c = 0, d = 0. 
Case 1.1.1: e = 0. 
In this case we have 

^112 ~ 0) ^122 ~ 0) ^222 ~ 0) ^111 ~ 0) ^112 ~ 0) ^122 ~ 0) (4-3) 

Al,, = -2a^ Al^^ = -2/2. (4.4) 

Hence all the A*^^ are zero except A[^^ and ^222^ which must be non-positive since a and / are 
real. 

Case 1.1.2: e 7^ 0. 

Simple manipulations yield 

^112 — 0) ^122 — 0) ^222 — 0) ^111 — 0) (4-5) 

= -2a\ Al,^ = -6a^ A\^^ = -6a/, Al^^ - -2j\ (4.6) 



which yield 



^ — ^111/2, /^ — ^222/^) 



{A\^^f - 9Aln^222 = 0, - SAl,, = 0. (4.7) 
This further requires ^in, ^n2; ^222 ^ 0- 

Case i.^: c = 0, d 0. 
Straightforward calculations lead to 

e'^-df-ae^ 0, (4.8) 

and 

= 0, ^^22 = 0, A^22 = 0, (4.9) 

= -2a^ A^ii = -2(a(i + de), (4.10) 

= -6e\ Al,^ = -6ef, Aj,^ = -2f\ (4.11) 
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which yield 



2 ' 2[-2>A\,,- Al,,±2^?>A[,,AU 

^ — ~^112/^J P — "^222/2, 



with conditions 



(^122) 3A]_;^2^222 — 0) 



^111 ^1 



^222 ^122 

Thus we require that AJii, A\^2 ^ 0; ^222 < 0- 

Case 1.3: c ^ 0, = 0. 
We find that e = a and 

^112 = 0) ^111 — 0) 
^111 — ~2a , ^122 — — 6ac, A222 — ~2c/, 
^^12 = -6a', A?22 = -6a/, >1^22 = -2(ac + f ), 



which yield 



2 _ ^111 2 _ (^122)^ f2 _ (^122 
" - O ' ^ - 1Q/I1 ' 



^112 "^^11 — 0) 9^111^222 3A];;^]_74j22 (^122) — 0) 



8l(>lin)'(^222)' - {A22)\Al,,f = 0. 

In this case we need to have A\^^ < 0, Al^2 ^ 0- 
Case 2: b^O. 

The subcases are as follows. 

Case 2.1: d = 0. 

In this case we find that 

e = 0, 6' - ac - 6/ = 0, 
9 



and 



^111 = 0, Al^^ = 0, Al^^ = 0, (4.22) 
A\^^ = -2a^ = -6a6, A}22 = -6&', (4.23) 
A^22 = -2(&c + c/), A222 = -2f, (4.24) 



which imply that 



= -Aiii/2, 6^ = -AI22/6, f = -^222/2, 

= ^ ^-!^^^^—^==^ (4.25) 

2[~3A|22 ~ ^122 ^ 2y3A];22^222] 



with extra conditions on A^^f^i given by 

{A\,^f - SAl,,A\,, = 0, (4.26) 



3AL2 = [±\/# + #][^l22 ± 73^122^122]- (4.27) 
" ^111 ^112 



and ^111 < 0, ^^22' ^222 — 



Case 2.2: d 0, e = cd/fe. 
We find that 

+ c^d - h^f - abc = 0, b^d + c^d^ - b^df - abed = 0, (4.28) 

and 

A\^^ = -2(a^ + bd), A\^^ = -6(a6 + be), A\^^ = -6(6^ + ce), (4.29) 

A^22 = -2(6c + c/), A^ii = -2(aci + cie), A^,^ = -6(60? + 6^), (4.30) 

^?22 = -6(te + e/), ^^22 = -2(ce + f ). (4.31) 

In order to construct the hnearizability conditions for above system one has to replace e = A| or 
e = A| to get a 5 X 5 matrix and the determinant of the resultant matrix must be zero. 
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5 Examples 



To test our linearization criteria, we utilize the following examples. 

1. For the case when at least one of the coefficients is variable e.g. P = 0, Q = 0, R = 3, S = 
0, T = 0, U = — = 0, W = 2, the system of two third-order equations is linearizable as 
it satisfies conditions (3.13)-(3.15) with a,b,c,d,e and / given in the Appendix. 



The following set of six equations can be obtained from [TO] by using the values of a, b, c, d, e and 
/ 

Px = 0, q^ = q/x, = 2r/x, 

Py = 2q/x, qy = —xp + r/x, Vy = —2xq. (5.1) 
A solution of system (5.1) is 

p = 1, q = 0, r = x^ . 
To obtain the Cartesian coordinates u{x,y) and v{x,y), one has to solve (JiO\) 

ul + vl = 1, 



UxUy + V^Vy = 0, 

ul + vl = x\ (5.2) 

Thus 

u = X COS y, 
t> = X sin y, 

are the transformations that linearize the system. 

2. The system which has P = 0, Q = -3, P = 0, ^ = 0, T = 0, U = 9y'^, V = 18, W = -6/y'^, 
is linearizable since it holds conditions (3.13)-(3.15) with a, b, c, d, e and / given in the 
Appendix. In the same manner as in example 1, one can obtain 

p = 1 + x'^ — 2x/y + l/y'^, 
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q = {l+x^)/y^ -x/y^, 



r = (l + x^)/y\ 

For the coordinates transformations one can solve system (5.2) with the help of p, q, r to obtain 

u = X- 1/y, 

V = x^/2 — x/y. 

3. The system of two third order equations with constant coefficients P = —2, Q = 0, R = 
—6, S = 0, T = 0, U = —6, = 0, = — 2 is also linearizable as the determinant of ( 14. ip is 
zero. 

The coordinate transformation for the above system which gives the linearization is 

u = cie-y-^ + C2ey-\ (5.3) 



V = cie"^-^ - cse^-^ (5.4) 

Note that there are four arbitrary constants which will appear in the solution for each of the above 
examples and therefore we have conditional linearizability subject to a system of two second-order 
equations. 



6 Concluding Remarks 



In this paper we have provided conditional linearizability criteria for a class of third-order systems 
of ODEs. To the best of our knowledge the linearization of systems of third-order ODEs has not 
been studied before in the literature. The system of two third-order ODEs that is cubically 
semi-linear in the first derivative was obtained by differentiating a system of two quadratically 
semi-linear ODEs and then we used the original system to replace the second derivative. The 
resultant system was then not a total derivative. 

The criteria developed are discussed in two parts. In the ffist part we constructed the linearizabil- 
ity conditions when at least one of the coefficients is a non-constant function, so that the system 
obtained for the linearizability conditions be invertible. For the constant coefficients case we 
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obtain a system which is independent of the derivatives. The non-trivial solution of the resulting 
system requires that the determinant of the coefficients be zero. There were various sub-cases. 
The linearizability conditions for each of the sub-cases of constant coefficients were constructed. 
This procedure will hopefully give rise to further studies in the construction of linearizability 
criteria for systems of higher-order equations. 
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Appendix 

By using Mathematica one can find a, b, c, d, e and / from system (3.15) in terms of 
P, Q, R, S, T, U, V and W as 

a = ^[3{-9S'^TU^ + RSU^ + 27S^TW - QSUW - SQSTV^ + SPSUV'^ 

-27RST^W + ISQSTUW - QRU^W + 3QRTVW - 27PSTVW + Q^UVW 
-3PRUVW - 9Q^TW^ + 21PRTW^){3Py - Q^)) + 2{9RSTU^ - R^U^ 
-27RST^V + 2QRU^V - 9PSU^V + 27PSTV'^ - Q^UV^ + 27R^T^W 
-18QRTUW + 9PRU^W + 9Q^TVW - 27 PRTVW){Qy - R^) 
+3{QT - PU){-3SU'^ + 9STV + RUV - QV'^ - 9RTW + 3QUW){Ry - 35^) 
-3{{QRSU^ - 9PS^U^ - 3QRSTV + 27PS^TV - Q^SUV + 3PRSUV 
+3QR^TW - 27PRSTW - Q^RUW - 3PR^UW + 18PQSUW + Q^VW 
-27 P^SVW - 9PQ^W^ + 27 P^RW^){3Ty - U^) + 2{QR^U'^ 
-9PRSU'^ - 9Q'^STV + 27 PRSTV - 2Q'^RUV + ld>PQSUV + gV^ 
-27 P'^SV'^ + 9Q'^RTW - 27 PR^TW - 9PQ^VW + 27 P'^RVW){Uy - K) 
-3{QT - PU){R^U - 3QSU - QRV + 9PSV + 3Q^W - 9PRW){Vy - 3W^)], (A.l) 

b - ^[-3{3S'^U^ - 9S'^TUV - 2RSU^V + 3RSTV'^ + 2QSUV'^ - 3PSV^ + 9RSTUW 
+R^U^W - GQSU^W - 3R^TVW + 9PSUVW - Q^V^W + 3PRV^W + 3Q'^UW^ 
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-9PRUW'^){3Py - Q^) - 2{RU - QV){-?>SU^ + 9STV + RUV - QV^ - 9RTW 
+3QUW)(Qy - R^) - (RU - QV){RU^ - 3RTV - QUV + 3PV^ + 9QTW 
-9PUW){Ry - 3S^) - 3{R'^SU^ - 3QS^U^ - 3R^STV + 9QS'^TV - Q^SV^ 
+3PRSV'^ + 3R^TW - 9QRSTW - 2QR^UW + GQ^SUW + 2Q^RVW - 3PR^VW 
-9PQSVW - 3Q^W'^ + 9PQRW^){3Ty - U^) + 2{RU - QV){R^U - 3QSU - QRV 
+9PSV + 3Q^W - 9PRW){Uy - 14) - {RU - QV){3R^T - 9QST - QRU 

+9PSU + g V - 3PRV) (Vy - 3W^)], (A.2) 

c = ^[-3{SV - RW){3SU'^ - 9STV - RUV + QV'^ + 9RTW - 3QUW){3Py - Q^) 

-6(-SV + RW)(RU^ - 3RTV - QUV + 3PV^ + 9QTW - 9PUW){Qy - R^) 
-{R^U'^V - 9QSTV'^ - 2QRUV^ + 9PSUV^ + QV^ - 9R^TUW + 27QSTUW 
-27 PSU^W + 18QRTVW - 9PQV^W - 27Q^TW^ + 27 PQUW^){Ry - 35^) 
+3{-SV + RW){R^U - 3QSU - QRV + 9PSV + 3Q'^W - 9PRW){3Ty - U^) 
-6(3R^T - 9QST - QRU + 9PSU + gV - 3PRV)(-SV + RW)(Uy - V^) 
+{-9R'^STU + 27QS^TU + R^U^ - 27PS^U^ - 2QR^UV + 18PRSUV + Q^RV^ 
-9PQSV^ + 9QR^TW - 27Q^STW - 9PR^UW + 27 PQSUW) {Vy - 3W^)], {A.3) 

d = ^[9RSTU'^ - R^U^ - 27 RST'^V + 2Qi?C/V - 9PSU'^V + 27 PSTV'^ - Q^UV'^ 

+27 R^T^W - 18QRTUW + 9PRU^W + 9Q^TVW - 27 PRTVW){3Py - Q^) 
+6{QT - PU){-3SU^ + 9STV + RUV - QV^ - 9RTW + 3QUW){Qy - R^) 
+3{QT - PU){RU^ - 3RTV - QUV + 3PV^ + 9QTW - 9PUW){Ry - 3S^) + QR'^U^ 
-9PRSU^ - 9Q^STV + 27PRSTV - 2Q^RUV + 18PQSUV + Q^V'^ - 27P^SV^ 
+9Q^RTW - 27 PR^TW - 9PQ^VW + 27 P^RVW){3Ty - U^) - Q{QT 
-PU){R^U - 3QSU - QRV + 9PSV + 3Q'^W - 9PRW){Uy - \4) 
+3{QT - PU){3R^T - 9QST - QRU + 9PSU + gV - 3PRV){Vy - 3W^)], {AA) 
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e = ^[-{RU - QV){-3SU'^ + 9STV + RUV - QV'^ - 9RTW + 3QUW){3Py - Q^) 

-2(RU - QV)(RU^ - 3RTV - QUV + 3PV^ + 9QTW - 9PUW){Qy - R^) 
-3(-R^TU^ + 3QSTU^ - 3PSU^ + 3R^T^V - 9QST^V + 9PSTUV 
+2PRU^V + Q^TV^ - 6PRTV^ - 2PQUV^ + 3P - 3Q^TUW 
+3PQU^W + 9PQTVW - 9P^UVW){Ry - 3S^) + {RU - QV){R^U - 3QSU 
-QRV + 9PSV + 3Q^W - 9PRW){3Ty - U^) - 2{RU - QV){3R^T - 9QST 
-QRU + 9PSU + gV - 3PRV){Uy - V^) + 3{3R^T^ - 9QRST^ - 2QR^TU 
+3Q'^STU + 9PRSTU + PR'^U'^ - 3PQSU'^ + 2Q'^RTV - QPR^TV - PQ'^V^ 

+3P'^RV'^ - 3Q^TW + 9PQRTW + 3PQ^UW - 9P^RUW){Vy - 3W^)], {A.5) 



f = —[S(SV- RW){RU'^ - 3RTV - QUV + 3PV'^ + 9QTW - 9PUW){3Py - Q^) 
-2(i?2f/V _ <^QSTV^ _ 2QRUV^ + 9PSUV'^ + Q^V^ - 9R^TUW + 27QSTUW 

-27 PSU^W + 18QRTVW - 9PQV'^W - 27Q'^TW'^ + 27PQUW^){Qy - R^) 
-3{-R^TUV + 3QSTUV - 3PSU^V + QRTV^ + PRUV^ - PQV^ + 9R^T^W 
-27QST^W - 3QRTUW + 27PSTUW + 3PRU^W - 18PRTVW + 9PV^W^ 
+27 PQTW'^ - 27 P^UW^){Ry - 3S^) - 3{3R^T - 9QST - QRU + 9PSU + QV 
-3PRV){-SV + RW){3Ty - U^) + 2{-9R^STU + 27QS^TU + R^U^ - 27PS^U^ 
-2QR^UV + 18PRSUV + Q'^RV^ - 9PQSV^ + 9QR^TW - 27Q^STW - 9PR^UW 
+27 PQSUW){Uy - V^) - 3(-9R^ST^ + 27QS^T^ + R^TU - 27 PS^TU - QR^TV 
-3Q'^STV + 18PRSTV - PR^UV + 3PQSUV + PQRV^ - 9P^SV^ + 3Q'^RTW 

-27PQSTW - 3PQRUW + 27 P'^SUW) {Vy - 3W^)] , (A.6) 

where 

A = (2{-R^{U^ - 27T^W) - 9Q\2ST + PW){V^ - 3UW) + 27 PS{SU^ - 3STUV 

-PV^ + 3PUVW) + Q=^(\/^ - 27TW'^) + 3R^{3ST{U^ - 3TV) + 6P{U^ - 3TV)W 

+QU{UV - 9TW) + 27QS{ST{-U^ + 3TV) + P{UV^ - 2U^W - 3TVW)) 
+3R{Q^V{-UV + 9TW) - 3Q{ST - PW){-UV + 9TW) + 9P{PW{V^ - 3UW) 

+S{-U^V + 2TV^ + 3TUW))))). (A 7) 
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